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Re normalization of Anisotropy and Glueball Masses on Tadpole Improved Lattice 

Gauge Action 

Mushtaq LoanQ Tim Byrnes, and Chris Hamer 
School of Physics, The University of New South Wales, Sydney, NSW 2052, Australia 

(Dated: January 10, 2003) 

The Numerical calculations for tadpole-improved 1/(1) lattice gauge theory in three-dimensions 
on anisotropic lattices have been performed using standard path integral Monte Carlo techniques. 
Using average plaquette tadpole renormalization scheme, simulations were done with temporal lat- 
tice spacings much smaller than the spatial ones and results were obtained for the string tension, 
the renormalized anisotropy and scalar glueball masses. We find, by comparing the 'regular' and 
'sideways' potentials, that tadpole improvement results in very little renormalization of the bare 
anisotropy and reduces the discretization errors in the static quark potential and in the glueball 
masses. 
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I. INTRODUCTION 

Compact U(l) gauge theory in (2+1) dimensions is one 
of the simplest models with dynamical gauge degrees of 
freedom and possesses some important similarities with 
QCD y|. The model has two essential features in com- 
mon with QCD, confinement |E!3 and chiral symmetry 
breaking 0. The theory is interesting in its own right, 
for it has analytically been shown to confine electrically 
charged particles even in the weak-coupling regime (at 
zero temperature) 0, IE IE IE • The confinement is 
understood as a result of the dynamics of the monopoles 
which emerge due to the compactness of the gauge field. 
The string tension as a function of the coupling behaves 
in a similar fashion to that of the 4-dimensional SU (N) 
lattice gauge theory. This model also allows us to work 
with large lattices with reasonable statistics. Other com- 
mon features of compact U (l)(2+i) and QCD are the exis- 
tence of a mass gap and of a confinement-deconfinement 
phase transition at some non-zero temperature. Thus, 
being reasonably simple and theoretically well under- 
stood in the weak-coupling limit, the U(l) model pro- 
vides a good testing ground for the development of new 
methods and new algorithmic approaches. In a recent 
paper we have obtained the first clear picture of the 
static quark potential, showing very clear evidence of the 
linear confining behaviour at large distances. The evi- 
dence of the the scaling behaviour of the string tension 
and the mass gap has also been observed in this model. 

In thepresent paper we want to extend the analysis 
of Ref. 9] in various respects. Since the measured ratio 
of spatial to temporal lattice spacings is not the same as 
the input parameter in the action, it becomes important 
to determine the true or renormalized anisotropy £ p hys. 
as a function of the bare anisotropy £o- An important 
advantage of using anisotropic lattices has been the need 
to measure the renormalization of anisotropy in the sim- 



ulation. The existing theoretical 0, and numerical 
studies E2 with Wilson action for SU(3) lattice gauge 
theory have shown that at finite coupling g, the renor- 
malized anisotropy £ p hys. differs appreciably form bare 
anisotropy £o. However, it has recently been shown that 
the use of improved actions, supplemented by tadpole 
improvement, besides providing a better discretization 
scheme for QCD, offer the advantage of a significant re- 
duction of renormalization of £o to a few percent |l3LIT^| . 
The anisotropy parameter r\ (ratio of the renormalized 
and bare anisotropy) and anisotropic coefficients have 
been calculated to one-loop order for improved actions in 
various recent studies [ll|, [IE HE HE 03 ■ These calcula- 
tions have provided very reliable results and the observed 
behaviour is confirmed non-pert urbati vely by large scale 
simulations on fine lattices El El HI We investi- 
gate the influence of tadpole improvement on isotropic 
and anisotropic lattices for the U(l) model in (2+1) 
dimensions in reducing the renormalization of the bare 
anisotropy at weak and strong couplings. We also apply 
tadpole improved U(l) lattice gauge theory to calcula- 
tions of the static quark potential, the string tension and 
the scalar glueball masses and compare the results with 
simulations of the Wilson action. 

The rest of this paper is organized as follows. Af- 
ter outlining the tadpole improved U(l) gauge model in 
(2+1) dimensions in Sect. II, we describe the method for 
determination of the renormalized anisotropy in Sect. III. 
We present results from our simulations on anisotropic 
lattices using both standard and tadpole improved Wil- 
son gauge action in Sect. IV. We compare bare and renor- 
malized anisotropies, static quark potential and scalar 
glueball masses from these actions. We conclude in Sect. 
V with a summary and outlook on future work. 



II. COMPACT 17(1) MODEL IN (2+1) 
DIMENSIONS 
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The tadpole-improved U{1) gauge action on an 
anisotropic lattice can be written in the following form 
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where P^ v is the plaquette operator, £o = At = 



(1) 

a t /a s 

is the bare anisotropy at the classical level and u s and 
Ut are the mean fields for the tadpole improvement. The 
notation used in Eq.Q differ slightly from that used in 
Refs. 113 ■ where the spatial and temporal mean- 
field improvement factors, u s and u t were absorbed into 
definition of [3 and £o- This, however, follows the notation 
introduced in Ref. Q. 

On the anisotropic lattice, the mean fields are deter- 
mined using the measured values of the average plaque- 
ttes [2l]]. We first compute u s from spatial plaquettes, 
u A s = (Pij), and then we compute u t from temporal pla- 
quettes, u1u 2 s = (Pit). Another way to determine the 
mean fields is to use the mean links in Landau gauge [22 



u t = (ReTrC/ t 



u s = (ReTrt/,; 



(2) 



E 



where the lattice version of the gauge condition is ob- 
tained by maximizing the quantity, 

^TTrlV(r) (3) 

Since the temporal lattice spacing in our simulations is 
very small, we adopt the following convention [L^ ITil. l2l| 
for the mean fields in tadpole improvement 

u t = l, Us = (P ZJ ) 1/4 . (4) 

This prescription eliminates the need for gauge fixing and 
the results yield values for u s which differ from those 
using Landau gauge by only a few percent. 



III. RENORMALIZATION OF ANISOTROPY 

Following the procedure of Klassen and Shake- 
speare and Trottier |19| . we measure the static quark 
potential extracted from Wilson loops in the spatial and 
temporal directions. Accordingly on an anisotropic lat- 
tice there are two potentials, V xt {R) and V xy (R). The 
two potentials differ by a factor of S, p hys. and by an ad- 
ditive constant, since the self-energy corrections to the 
static potential are different if the quark and anti-quark 
propagate along the temporal or a spatial direction. Thus 
£phys. can be determined by comparing the static quark 
potential computed from the logarithmic ratio of time- 
like Wilson loops R(x,t), where 



R(x,t) 



_ W xt (x,T + l) 



(5) 



W xt { X ,T) ' 

with the potential computed from that of the space-like 
Wilson loops R(x,y), where 

W xy (x,y + 1) 



R{x,y) 



W xt (x,y) 



(6) 



Asymptotically, for large r and y, the ratios R(x,t) and 
R(x,y) approach 



R(x,t) = Z XT c 

>:>J 



-rV xt 



(excited state contr.) (7) 
R{x,y) = Z xy c~ yV:cy + (excited state contr.). (8) 



To suppress the excited state contributions, a simple 
APE smearing technique [2^, |3 |2f| was used. In this 
technique an iterative smearing procedure is used to con- 
struct Wilson loop (and glueball) operators with a very 
high degree of overlap with the lowest-lying state. In our 
single-link smoothing procedure, we replace every space- 
like link variable by 



U, 
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(9) 



where the sum over V refers to the "staples", or 3-link 
paths bracketing the given link on either side in the spa- 
tial plane, and P denotes a projection onto the group 
C/(l), achieved by renormalizing the magnitude to unity. 
We used a smearing parameter a = 0.7 and up to ten 
iterations of the smearing process. To reduce the statis- 
tical errors, the time-like Wilson loops were constructed 
from "thermally averaged" time- like links [2J, |2!| . 

The links making up the space-like and time-like Wil- 
son loops are smeared by the same amount so that the 
ratios R(x,t) and R(x,y) have the same excited-state 
contribution. Similarly, finite-volume corrections to the 
R(x, t) and R(x, y) are the same if the temporal and spa- 
tial extents are equal in physical units, i.e. N s = ^phys.Nt 
in lattice units. These statements are expected to hold 
only for large x, y and t; otherwise there can be large 
0(al,af) lattice errors. 

The physical anisotropy is determined from the ratio 
of the potentials V x t (R) and V xy (R) estimated from R xt 
and R xy respectively. The unphysical constant in the 
potentials is removed by subtraction of the simulation 
results at two different radii 



£,phys 



V x t{R2)-V xt (Ri) 

V X y(R2) — V X y(R\) 



(10) 



The measured renormalization of the anisotropy, 77 is then 
determined from 



Cphys 



(11) 



IV. SIMULATION AND RESULTS 



Simulations were performed on four lattices of N 2 x Nt 
sites, with N s — 16 and Nt ranging from 32 to 48 with 
mean-link improvement, and four lattices with the Wil- 
son action. Configurations were generated by using the 
Metropolis algorithm. The details of the algorithm are 
discussed elsewhere @ . 50000 sweeps were performed for 
thermalization of the configurations and self-consistent 
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determination of the tadpole factors. Configurations are 
stored every 250 sweeps thereafter. Ensembles of about 
1000 configurations were used to measure the static quark 
potential, while 1,400 configurations, at coupling values 
from (3 — 1.0 to 2.5, were generated for the glueball mass. 
We fixed £o = 16/iVt in the first pass, so that the lattice 
size remains fixed at 16a s in all directions. The simula- 
tion parameters of the lattices analyzed here are given in 
Table D 



TABLE I: Simulation parameters for the lattices analyzed 
for renormalization of the anisotropy. The bare anisotropies 
£o and the means fields u t and u s for tadpole improvement 
are shown. 



Action 


Co = Ar 


P 


u t 




Volume 


Tadpole 


0.50 


1.306 


1. 


0.924 


16 2 x 32 


Improved 




1.443 


1. 


0.931 








1.592 


1. 


0.937 








1.892 


1. 


0.946 






0.40 


1.302 


1. 


0.921 


16 2 x 40 






1.436 


1. 


0.927 








1.584 


1. 


0.932 








1.88 


1. 


0.940 






0.33 


1.2931 


1. 


0.915 


16 2 x 48 






1.426 


1. 


0.920 








1.567 


1. 


0.922 








1.858 


1. 


0.929 




Wilson 


0.50 


1.4142 






16 2 x 32 




0.40 








16 2 x 40 




0.33 








16 2 x 48 



0.2 



AT = 0.5 



AAA -A— — — ^ — X~ 



FIG. 1: Logarithmic ratio of the time-like Wilson loops as 
a function of r for fixed x at (5 = 1.309 and At = 0.5. The 
solid triangles correspond to x = 4 and open triangles show 
x = 2 values. 
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After measuring the Wilson loops at fixed values of (3, 
we compute the ratios R x t and R xy . We find that the 
individual ratios reach their plateaus for r > 3 and y > 3 
for fixed x as shown in Figures ^ and [21 These ratios are 
are expected to be independent of r and y for r, y > 3 
respectively. The estimates of the potentials V x t (R) and 
V xy {R) can now be found from these ratios. 

Figure shows a graph of the static quark potentials, 
computed from spatial and temporal Wilson loops, as a 
function of radius R at f3 = 1.306 and At — 0.5. The 
potential in the lattice units obtained from the ratio of 
the time-like Wilson loops has been rescaled by the input 
anisotropy. To extract the string tension, the time-like 
potential is well fitted by a form 

V xt (R) =a + bliiR + aR, (12) 

including a logarithmic Coulomb term as expected for 
classical QED in (2+1) dimensions which dominates the 
behaviour at small distances, and a linear term as pre- 
dicted by Polyakov Q and Gopfert and Mack 8] dom- 
inating the behaviour at large distances and showing a 
clear evidence of the linear confining behaviour at large 
distances. 



0.25 1 1 1 1 1 1 1 1 1 

012345678 

y 

FIG. 2: Logarithmic ratio of the spatial Wilson loops as a 
function of y for fixed x at /? = 1.309 and At — 0.5. The 
solid triangles correspond to x — 4 and open triangles show 
x — 2 values. 



We measured each anisotropy twice, using two different 
radii R\ for subtraction, with fixed i?2- Setting R2 = 
4, we computed the anisotropy with R\ = 2 and 
The two determinations of anisotropy, shown in Table ITT1 
are in excellent agreement. The numerical values of the 
renormalization of the anisotropy parameter 77 appears 
to be equal to unity even at large f3. ft is seen that 
for with mean-field improvement the input anisotropy is 
renormalized by few percent over the range of lattices 
analyzed here, whereas the measured value of anisotropy 
is about 15 — 20% lower than the bare anisotropy with the 
standard Wilson action. This can be seen from Figure 
2] where the renormalization of the anisotropy is plainly 
visible as a difference in slope of the potentials computed 
from R xt and R xy . 

Glueball correlation functions C(r) were also calcu- 
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V(R) 




FIG. 3: Static quark potentials computed from R x t (solid 
triangles) and R xy (open squares) as a function of separation 
R at P — 1.306 and At = 0.5 using mean field improved gauge 
action. The potential obtained from R xt has been rescaled 
by the input anisotropy. The solid line is a fit to the form 
V(R) = a + ctR + c\n(R) to the temporal Wilson loops. 



V(E) 



FIG. 4: Static quark potentials computed from R x t (solid 
triangles) and R xy (open squares) at P = 1.306 and At — 0.5 
for standard Wilson action. 



lated 



C(r) = ^(0|^(r + To)^(ro)|0) 



(13) 



where $j(r) is the optimized glueball operator found by 
a variational technique, following Morningstar and Pear- 
don and Teper [2(| , from a linear combination of the 
basic operators (pi, 



$(r) 



Via4>ia{r) 



(14) 



where the index a runs over the rectangular Wilson loops 



with dimensions 
and smearing n s — 
in all. 



= [n — l,n + 1], l y = [n — 1, n + 1] 
[m — 1, m + 1], making 27 operators 



TABLE II: Measured anisotropy £. p h y s. compared to the bare 
anisotropy £o for two actions. 



Action £q p 



Ri 



£,phys. Tj — £,phys./^0 

2 iii = V2 Ri = 2 Ri = V2 



Tadpole 0.50 1.306 0.500(1) 0.500(2) 

Improved 1.443 0.501(2) 0.503(2) 

1.592 0.50(2) 0.50(3) 

1.892 0.50(5) 0.50(4) 

0.40 1.302 0.402(2) 0.401(2) 

1.436 0.402(1) 0.403(8) 

1.584 0.40(3) 0.40(2) 

1.88 0.40(1) 0.40(1) 

0.33 1.293 0.335(3) 0.34(1) 

1.426 0.33(2) 0.33(1) 

1.567 0.33(2) 0.33(1) 

1.858 0.33(3) 0.33(1) 

Wilson 0.50 1.414 0.482(4) 0.471(3) 

0.40 1.414 0.336(4) 0.338(3) 

0.33 1.414 0.279(5) 0.281(2) 



1.00(1) 
1.00(4) 
1.00(4) 
1.0(1) 
1.00(4) 
1.00(2) 
1.00(6) 
1.00(7) 
1.00(8) 
1.00(6) 
1.00(6) 
1.00(6) 



1.00(4) 
1.00(4) 
1.00(6) 
1.00(8) 
1.00(5) 
1.00(1) 
1.00(4) 
1.00(2) 
1.03(3) 
1.00(3) 
1.00(3) 
1.00(3) 



0.964(8) 0.942(6) 
0.84(1) 0.845(7) 
0.83(1) 0.843(6) 



The optimized correlation function was fitted with the 
simple form 



Ci — c + ci cosh mi(T/2 — r) 



(15) 



to determine the glueball mass estimates. 

The results for the symmetric and the anti-symmetric 
glueball masses over the square root of the string tension 
are shown in Table II I II along with the mean plaquette 
values at different (3 at At = 1.0. Figure shows the be- 
haviour of the logarithm of the antisymmetric mass gap 
over the square root of the string tension as a function of 
(3. It can be seen that that the ratio scales exponentially 
to zero in the weak-coupling limit as should be in three- 
dimensional confining theories. The solid line is a fit to 
the data over the range 0.916 < (3 < 2.12. The slope 
of the data matches the predicted form pj, however, the 
intercept of the scaling curve is large by a factor of 2 
(our previous estimates of constant coefficient for stan- 
dard Wilson action are large by a factor of 5.2 0). It 
would be interesting to test the sensitivity of the slope 
and intercept of the scaling curve by including the radia- 
tive corrections. 

A plot of mass ratio against the effective lattice spacing 
a eff is shown in Figure El At weak coupling, the 
theory is expected to approach a theory of free bosons 
Q so that symmetric state will be composed of two 0~" 
bosons and the mass ration should approach two. The 
mass ratio of the lowest glueball states scale against the 
effective lattice spacing towards a value close to 2.0, as 
expected for a theory of free scalar bosons. We note 
that mass ratio exhibits scaling behaviour, even with the 
Wilson action however, in contrast with the Wilson 
action, a significant reduction in the errors with tadpole 
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cq. 1 




1.2 1.4 



P 



FIG. 5: Logarithmic dimensionless ratio M/V^/3 3//4 as a 
function of f3. The solid curve is the fit to the data for 0.916 < 
/3 < 2.12. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 6: Mass ratio Mr against the effective lattice spacing 
a e ff at At = 1.0. The solid triangles show the estimates for 
mean field improved gauge action and open circles show the 
estimates of standard Wilson action. 



improved action is apparent in the mass gap and the mass 
ratio. 



V. SUMMARY AND OUTLOOK 

Mean field improved U{\) lattice gauge theory in (2+1) 
dimensions was applied to calculations of the static quark 
potential, the renormalized lattice anisotropy and the 
scalar glueball masses. We analyzed the mean-link im- 
proved action on isotropic and anisotropic lattices and 
comparisons were made with the simulations of the Wil- 
son action. By comparing the static quark potentials 
computed from space-like and time-like Wilson loops, we 
determined the physical anisotropy of the tadpole im- 
proved Wilson action. We found that with mean-link 
improved Wilson action, the bare anisotropy is renor- 



malized by less than a few percent, in contrast with the 

TABLE III: Monte Carlo estimates for the mean plaquette 
(P) , the effective lattice spacing, the symmetric and antisym- 
metric glueball masses M ++, M --, for tadpole improved 
Wilson action at At = 1.0 



Action 


13 


<^> 


a eff 


m ++/Vk 


m —/Vk 


Rm 


Tadpole 


0.917 


0.678 


0.864 


4.0(1) 


3.9(1) 


1.03(3) 


improved 1.243 


0.719 


0.446 


3.83(9) 


3.2(1) 


1.19(3) 




1.327 


0.776 


0.373 


3.5(2) 


2.9(1) 


1.23(4) 




1.460 


0.797 


0.281 


3.7(1) 


2.34(7) 


1.58(4) 




1.613 


0.816 


0.202 


3.6(1) 


2.05(2) 


1.80(5) 




1.816 


0.836 


0.129 


2.05(4) 


0.98(3) 


2.00(5) 




1.917 


0.844 


0.103 


2.0(1) 


0.99(5) 


2.00(6) 




2.168 


0.861 


0.058 


2.19(6) 


1.05(2) 


2.08(7) 




2.417 


0.874 












2.669 


0.885 












2.917 


0.895 










Wilson 


1.0 


0.475 


0.733 


3.27(3) 


3.1(1) 


1.03(6) 




1.35 


0.629 


0.356 


3.6(1) 


2.81(9) 


1.30(5) 




1.41 


0.656 


0.310 


3.8(1) 


2.95(8) 


1.31(7) 




1.55 


0.704 


0.231 


4.0(1) 


2.81(6) 


1.42(7) 




1.70 


0.748 


0.166 


4.03(8) 


2.51(2) 


1.60(8) 




1.90 


0.790 


0.107 


3.9(1) 


1.88(3) 


2.0(1) 




2.0 


0.806 


0.085 


3.53(4) 


1.96(4) 


1.80(9) 




2.5 


0.854 


0.027 


3.1(2) 


1.50(8) 


2.0(1) 



standard Wilson action, where the measured value of 
anisotropy is found to be about 15 — 20% lower than bare 
anisotropy on the lattices analyzed here. We found that 
tadpole improvement significantly reduces discretization 
errors in the static quark potential and the glueball 
masses. The mass ratio of the two lowest glueball states 
scales against the effective lattice spacing towards a value 
close to 2.0, as expected for a theory of free scalar bosons. 

We intend to extend PIMC techniques to Symanzik 
improved U(l) lattice gauge theory. The intention is to 
study the effects of improvement on the scaling slope, 
the constant coefficients and scaling behaviour observed 
in the weak-coupling regime of the theory. We also plan 
the study the one-loop correction to the anisotropy fac- 
tor for Symanzik improved U{1) gauge action in three 
dimensions. We shall report on this work in the near 
future. 
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